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1. Introduction 



Obtaining equilibrium properties for the hydrogen atom in free space, through standard 
Boltzmann-Gibbs (BG) statistical mechanics, is troublesome, since the partition function di- 
verges for any finite temperature. This occurs mostly in systems of composite particles, which 
are characterized by upper bounds, preceded by a quasi-continuum of energy levels, in their 
energy spectra. Each of the levels in the quasi-continuum yields a small contribution for the 
partition function; the divergence arises since there are, in principle, an infinite number of such 
levels. Essentially, this reflects the fact that composite particles always ionize at any finite 
temperature. Due to these difficulties, such systems are never discussed in standard statistical- 
mechanics textbooks (for an exception of this, see Ref. Q]). 

However, the presence of long-living nonionized hydrogen atoms in galaxy peripheries and in 
intergalactic media, at low temperatures, is undeniable. Certainly, their density must be very 
low in such a media, in order to avoid combination (H + H — > H2), and they should spend a long 
time in their nonionized state, before reaching ionization. These long-living quasi-equlibrium 
states represent the main point we explore, on theoretical grounds, in the present work. 

The specific heat of the hydrogen atom has been calculated recently j^l through the gener- 
alized statistical-mechanics formalism that emerged from Tsallis's generalization of the 
BG entropy Within this approach, the specific heat was computed for certain values of the 
entropic index q < 1 (notice that q = 1 corresponds to the standard BG formalism). Even within 
such a formalism, the hydrogen-atom specific heat presents several anomalies, like divergences, 
cusps, and discontinuities in its derivative P). 

Herein, we choose a different approach for analyzing this problem, as described next. 

(i) We investigate the dynamical behavior of the hydrogen atom by applying a standard Monte 
Carlo method, in which the probability for jumping between states is based on the BG weight. As 
expected, for any finite temperature, the simulation always carries the system towards ionization, 
after some time. However, it is shown that the system may live in a quasi-equilibrium state, 
characterized by a slowly varying value of the average energy, before approaching its maximum- 
energy state. The time that the system remains on such a quasi-equilibrium state increases for 
lowering temperatures. 

(ii) We propose a modified dynamics that prevents the system from reaching ionization and 
whose results coincide, during some time (essentially when the system lies in its low-temperature 
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quasi-equilibrium state), with those of the standard dynamics. Therefore, for low temperatures, 
depending on time scale of interest, the modified dynamics may reflect the correct dynamical 
behavior of the system. The advantage of the modified dynamics is that its corresponding 
partition function, associated to the statistical weight that generated it, is finite and may be 
calculated exactly. 

(iii) Since the quasi-equilibrium state may present a long duration for low temperatures, 
herein we will consider it as an effective equilibrium. Therefore, if the results obtained from the 
standard and modified dynamics coincide within such a quasi-equilibrium state, one may use 
the later formalism in such a way to compute thermodynamic properties. 

(iv) We introduce a modified regularized partition function, by extracting the divergence of 
the BG partition function, for a hydrogen atom, and compare the corresponding internal energy 
and specific heat with the results obtained from the numerical simulations. 

(v) Obviously, the procedures described above may work as good approximations for low 
energies, but should fail for increasing energies. Since our energies are always measured with 
respect to the corresponding ionization energy, we show that our approximations work well for 
energies that correspond to temperatures much higher than room temperature. 

In the next section we investigate the dynamical behavior of the hydrogen atom within a 
standard Monte Carlo framework. In section 3 we introduce a modified dynamics and apply 
it for the hydrogen atom. In section 4 we calculate a regularized partition function, related to 
the modified dynamics introduced in section 3, and compute some associated thermodynamic 
quantities for the hydrogen atom. In section 5 we propose two possible physical realizations for 
the quasi-equilibrium states presented herein. Finally, in section 6 we present our conclusions. 

2. The Hydrogen Atom within a Standard Monte Carlo Procedure 

Let us consider a hydrogen atom with its well-known energy spectrum (see, e.g., ref. [3], 



where R is the Rydberg constant [R = 13.6058 eV, or R = 2.18 x 10 18 J, which corresponds 
to a temperature {R/ks) = 1.579 x 10 5 K], and we have chosen the ground state to have zero 
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energy. As n increases, one has a quasi-continuum of energy levels, and in the limit n — * oo the 
atom ionizes (with an ionization energy E* = R), in such a way that the gap separating the 
ground-state and ionization energies is R. A transition between the ground state and the first- 
excited state costs most of the energy of the gap, i.e., 3R/4, whereas all further jumps occur in 
the range R/A. Therefore, for low temperatures, the hydrogen atom is expected to remain a long 
time in the lowest-energy states (mostly in the ground state). However, as the quantum number 
n increases, the energy cost for jumps between nearest-neighbor levels decreases. Therefore, 
once the system has reached a state characterized by a large quantum number n, transitions to 
higher-energy levels cost very little energy, in such a way that after a long time, the hydrogen 
atom will ionize. The ionization will always occur for any finite temperature; one of the questions 
we address in the present work is how long does the system remains in the quasi-equilibrium 
state (characterized by an average energy close to the ground-state energy). 

It is important to remind that the above energy spectrum is degenerate (see, e.g., ref. Q], 
page 798), in such a way that a given energy level E n presents a total degeneracy g n = n 2 . 
Therefore, a precise analysis of this problem should take into account the degeneracy of the 
energy levels, as well as of the spins. However, in the present investigation, only the lowest 
energy levels will contribute most significantly, in such a way that the effects of the degeneracy 
will not change qualitatively our results. Consequently, we will consider herein, g n = 1 (Vn), for 
simplicity. 

We have investigated the dynamical behavior of the hydrogen atom through a Monte Carlo 
procedure 8]. The probabilities P{n + 1 — > n) and P(n — ► n + 1), for transitions between the 
states characterized by quantum numbers n and n + 1, satisfying the detailed-balance condition 
and constructed by using the standard BG weight, are given by 

P(n + 1 -» n) = c , P(n -> n + 1) = c exp[-0(E n+1 - E n )\ , (2) 

where (3 = (fcgT)" 1 and c is an arbitrary constant (0 < c < 1/2). In order to implement the 
dynamical evolution, a uniform random number z (z £ [0, 1]) must be generated at each Monte 
Carlo step (which will be adopted as our unit of time). For a system in a state characterized by 
the quantum number n at time t, transitions between states are performed (or not), depending 
on the value of z(t), according to the following rules: 

(i) If z(t) < P(n — ► n + 1), perform the jump n — > n + 1; 

(ii) If P(n — > n + 1) < z(t) < [P(n — » n + 1) + P(n — * n — 1)], perform the jump n—*n—l] 
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FIG. 1: The time evolution of the average dimensionlcss energy (energy in units of the Rydberg constant), 
for the hydrogen atom, is represented for several relative temperatures (6 = ksT/R). Time is measured 
in Monte Carlo steps. 

(iii) Else, remain on level n. 
One may easily see that the constant c is proportional to the probability for no jumps [rule 
(hi)]. In the following results we have used c = 1/2, although we have verified that other choices 
for this constant (0 < c < 1/2) did not change qualitatively our results. Below, ( ) correspond 
to averages over N s distinct samples, i.e., different sequences of random numbers. We have 
considered two distinct initial conditions in our simulations, which will be reffered to, herein, 
as conditions 1 and 2, corresponding to all samples starting with the quantum number n = 1 
(condition 1), and all samples starting with n = 2 (condition 2). Although, in some cases, the 
results obtained by using these two initial conditions seem to be different, they approach each 
other in the limit N s — ► oo (except, of course, in the transient regimes, before reaching the 
quasi-equilibrium state) . 

In Fig. 1 we present the time evolution of the average dimensionless energy, (e(t)) [e n = 
E n /R], for several values of the relative-temperature variable, = (ksT)/R, i.e., the ratio of 
the temperature with respect to the Rydberg constant. Our simulations were carried up to a 
maximum time i max = 10 9 , whereas for the averages we have considered N s = 2000 samples 
with the initial condition 2. Our plots exhibit a general tendency for increasing the average 
energy, towards the ionization energy, after some time. This reflects the fact that the hydrogen 
atom always ionize for any finite temperature. However, the interesting effect noticed herein 
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FIG. 2: The two-time autocorrelation function C(t+t w , t w ), for the hydrogen atom, is represented against 
time, for a given value of the relative temperature and several waiting times t w . Time is measured in 
Monte Carlo steps. 

is the presence of a quasi-equilibrium state, characterized by a slowly varying average energy, 
before the approach to ionization. Such a quasi-equilibrium state may present a long duration 
for low values of 9, as shown in Fig. 1. 

It is important to notice that such a state presents quite nontrivial behavior, e.g. it is 
characterized by aging, as observed on a similar system |9|. In order to see this effect, let us 
define the two-time autocorrelation function 



EH 



r U-L* * ^ - ^ + ^) £ (*w)) - (e(* + t w ))(e(t w )) , , 

o"(* + t w )cj(t w ) 

where t w represents the well-known "waiting time" , and 

a(t) = [((e(t)f) - (e(t)) 2 ] 1/2 ■ (3b) 

In Fig. 2 we exhibit such a correlation function for a value of the relative temperature 9 = 0.08 
and several waiting times. The initial conditions are the same as in Fig. 1, but now, the averages 
that appear in Eqs. (3) were taken over N s = 350000 samples. The waiting times considered 
were chosen in such a way to ensure that the correlation function C(t+t w , t w ) was evaluated with 
the system on the quasi-equilibrium state (which is typically inside the time range 10 2 -> 10 8 , 
as shown in Fig. 1). Clearly, there is a dependence on the waiting time typical of the aging 
effect Q- 
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FIG. 3: The time evolution (in Monte Carlo steps) of the average dimensionless energy (energy in units of 
the Rydberg constant), for the hydrogen atom, is represented for two relative temperatures (9 = ksT/R), 
in the case of the standard (empty symbols) and modified dynamics (full symbols). 

3. The Hydrogen Atom and the Modified Monte Carlo Procedure 

We have also considered the evolution of the system under a modified dynamics, which 
satisfies detailed balance, but prevents the approach to ionization (the justification for that will 
become clear later on). The jumping probabilities are given by 

P(n+l-n) = c, P(n -> n + 1) = c " ^"ffi • (4) 

V ' V ; exp(-pE n ) - exp(-/3i?) v ; 

In Fig. 3 we compare the results for the time evolution of the average dimensionless energy 

obtained by using the standard dynamics [Eq. (2)] with those obtained through the modified 

dynamics [Eq. (4)], for typical values of the relative temperature. The conditions for the 

simulations in Fig. 3 are the same dbS those of Fig. 1, which correspond, to ^max — 

10 9 and 

N s = 2000 (initial condition 2). One observes that the agreement between the two dynamical 
procedures persists during a given time (essentially within the quasi-equilibrium state), that 
increases for lowering values of the temperature. Within the modified dynamics, which was 
constructed in such a way to avoid the system from reaching ionization, the system remains on 
a quasi-equilibrium state forever. 

Let us define the duration of the quasi-equilibrium state, £ms(#)> as the time during which 
the system remains on such a state, within the standard dynamics, by keeping the absolute value 
of the difference between the average dimensionless energies computed from the standard and 
modified dynamics less than a given value 5. Although the choice of 5 may be arbitrary, one 
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expects the corresponding law followed by <ms(^) to be independent of this particular choice. 
In the present analysis we estimated tyis(0) by considering several values of 6 from 0.08 up to 
0.12, by imposing that 5 does not exceed 0.05. Our data fit well the exponential law, 



*ms(0) ~ exp(6/0) = exp[bR/(k B T)] (b = 1.95 ± 0.06), (5) 
which implies tyis(Q) — > oo when T — > 0. Essentially, the duration of the quasi-equilibrium state 
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follows an Arrenhius law, typical of the Kramers' escape problem in chemical reactions 
where the system may remain for a long time in a quasi-stationary state, before overcoming the 
potential barrier associated with the reaction. To our knowledge, this is the first time that such 
a behavior has been associated with the dynamics of the hydrogen atom. 

Therefore, for very low temperatures and depending on time scale of interest, the modified 
dynamics may reflect the correct dynamical behavior of the system. The advantage of the 
modified approach is that the corresponding partition function, associated to the statistical 
weight that generated Eq. (4), is finite - contrary to the one associated with the standard 
dynamics - and may be calculated exactly. This will be done in the next section. 



4. The Regularized Partition Function and Associated Thermodynamic 
Functions 



In the present section we shall calculate the partition function defined through the statistical 
weight related to the modified dynamics, introduced above. Let us address this point by con- 
sidering in detail the divergence of the partition function associated with the energy spectrum 
ofEq. (1), 

oo oo 

Z = ^ exp(—/3E n ) = exp(- ^i?)^exp[W(n 2 )] . (6) 

n=l n=l 

The equation above may still be written as 

n* 

Z = exp(-/3i?) lim V exp[/3i?/(n 2 )] , (7) 

n*— >oo 

n=l 

where n* is an appropriated cutoff in the quantum number, that will be taken to infinite later 
on. One has that 
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(8) 



where H n * 2m = Xm=i l/( n2m ) are the harmonic numbers of order 1m |l3j |. The limits n* — > oo 
of the harmonic numbers are well-defined, leading to finite coefficients, Bim = lim n *_»oo Hn*,2m- 
One has that B 2 = vr 2 /6 = 1.64493..., B 4 = vr 4 /90 = 1.08232..., B 6 = vr 6 /945 = 1.01734..., 
in such a way that Bi m converges to unit for increasing values of m, e.g., B\q = 1.00002.... 
Therefore, one gets 



Z = exp(-0R) 



lim (n*) + V B 



2in 



(PRY 



m=l 



nil 



(9) 



which shows a linear divergence with the quantum number. It should be stressed that the 
divergent contribution of the partition function comes from a single term in the sum over m of 
Eq. (8) (term m = 0). Let us now introduce a "modified regularized partition function", 



Z' = Z- exp(-0R) lim (n*) = exp(-/?B) V B 



2m 



{my 



m=l 



ml 



(10) 



which is finite. The regularized partition function defined above may be written also as 
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Z' = J2^M-PE n ) - ^exp(-/3B) = [<sxp(-pE n ) - exp(-/3B)] , 

n=l n=l n=l 

where one identifies the statistical weight that leads to the jumping probabilities of Eq. (4). 

Obviously, Z and Z' are very distinct from one another (actually, the difference between 
them diverges). However, for low temperatures, the hydrogen atom remains for a long time in 
its low-energy states. Therefore, for a given low temperature, if one considers the corresponding 
quasi-equilibrium state as an effective equilibrium, for which the divergent term is not relevant, 
one may use Z' in order to calculate thermodynamic properties as approximations. Within this 
formalism, the quantities that would correspond to the internal energy and specific heat are 
given, respectively, by 
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FIG. 4: The internal energy of the hydrogen atom is calculated up to T = 15000 K, by two different ap- 
proximations: (i) The analytic modified regularized partition function formalism (full line); (ii) Standard 
Monte Carlo simulation (black dots). The inset shows an amplification of these results in the temperature 
range T = 0K^T = 11000 K. 
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where we keep the prime notations to remind that such thermodynamic properties were calcu- 
lated by using the regularized partition function Z'. 

In Fig. 4 we exhibit the internal energy of the hydrogen atom, as calculated from the above 
analytic expression [Eq. (12a)], and compare it with results of the numerical simulations. In the 
latter, we have computed time averages over a time interval ranging from t\ = 200 up to t 2 = 400 
Monte Carlo steps of the standard-dynamics quasi-equilibrium states, where the values at each 
time do already represent averages over N s = 10 6 samples, with the initial condition 1. Similar 
results may be obtained by starting the system with the initial condition 2, although a larger 
computational effort (i.e., higher values of N s ) may be necessary for a proper convergence towards 
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the low-temperature quasi-equilibrium state in this case. One observes a very good agreement 
between the two approaches up to temperatures 10 4 K. In fact, the relative discrepancy between 
these procedures gets larger for increasing temperatures, yielding the typical values of 0.004, for 
9 = 0.05 (T ps 7895 K) and 0.282, for 6 = 0.06 (T ps 9474 K). A similar picture to the one shown 
in Fig. 4 holds for the specific heat of the hydrogen atom, with relative discrepancies of the same 
order of magnitude as those found for the internal energy. For the range of temperatures over 
which one finds a good agreement between the two approaches (typically from K to 10 K), the 
long-living low-temperature quasi-equilibrium state may be considered as an effective equilibrium 
state, for which the divergence of the partition function in Eq. (9) has been removed, and one 
may compute effective thermodynamic properties for the hydrogen atom from the regularized 
partition function of Eq. (10). It is important to mention that the range of temperatures over 
which the present regularized formalism should be applied safely for the hydrogen atom goes far 
beyond room temperature. 



5. Possible Physical Realizations 



It seems difficult from the experimental point of view to perform measurements on a system 
of highly diluted nonionized hydrogen atoms, since one has to avoid the atoms from reaching 
ionization (which is favored at high temperatures), as well as from achieving combination, H + 
H — > H2 (which becomes enhanced at low temperatures). In the following discussion we propose 
two possible physical situations in which hydrogen atoms may be found in the above-mentioned 
long-living quasi-equilibrium state. In both cases, the duration of such a state is estimated in 
real time. 



A. Gas of Hydrogen Atoms 



Let us assume the viability, from the experimental point of view, for generating a gas of 
nonionized hydrogen atoms at low temperatures (as compared with the corresponding ionization 
temperature). We shall also assume that the following conditions are satisfied: 

(i) The combinations, H + H — * H2, are negligible. In fact, there are experimental techniques 
for such a purpose, in which pairs of atoms in the so-called "spin-polarized state", can not 



produce bound states 
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(ii) The temperature T is sufficiently low in such a way that most of the atoms that compose 
the gas are in the ground state; 

(iii) A given hydrogen atom can only change its state through collisions with other atoms 

(we will assume that the linear dimension of the box enclosing the gas is much larger than its 

mean free path, in such a way that collisions with the walls may be neglected). An atom in its 

ground state may experience several collisions before excitation (actually, such an atom should 

undergo a collision with a sufficiently energetic atom in such a way to absorb an energy greater 

than 3i?/4 for a transition to occur). 

Therefore, the mean time between two successive collisions, r, may be obtained from standard 

□ 

kinetic-theory calculations (see, e.g., ref. |17l j. chapter 12), 



1 / m \ 1/2 



where n represents the density of atoms in the gas, m is the mass of the hydrogen atom, and 
a stands for the Bohr's radius. In the present physical system, a given atom can only change 
its energy through collisions with other atoms; this situation is mimicked in the Monte Carlo 
simulation, where at each step there is a finite probability for the occurence of a given energetic 
transition. In order to establish a connection between the duration of the quasi-equilibrium 
states of Fig. 1 (given in Monte Carlo steps) with real time, we shall propose the crude - but 
very suggestive - correspondence: r = 1 Monte Carlo step. With this assumption, one gets that 
the duration (in real time) of a hydrogen atom quasi-equilibrium state, may be written as, 



1 / m \ 1/2 

ireal = *ms t « — -= £ TT~^r exp[2i?/ (ksT)] , (14) 



where we have considered the fitting parameter 6 = 2 [cf. Eq. (5)]. It should be stressed that, at 
low temperatures, the exponential growth dominates completely, with the multiplicative factors 
in Eq. (14) becoming irrelevant. Therefore, the hypothesis above for the connection of Monte 
Carlo steps with real time may be softened in such a way that any linear relationship between 
a given number of successive collisions and another number of Monte Carlo steps would lead to 
the same exponential-growth behavior of Eq. (14). 
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One sees that the above duration time depends on two parameters, namely, the density 
of atoms and the temperature, i.e., t rea i = t iea i(n,T). We have estimated t rea \(n, T) for the 
typical values n = 10 23 atoms/m 3 and 9 = 0.10 (T pa 15797 K), 9 = 0.05 (T Pa 7898 K), 
9 = 0.03 (T pa 4739 K), and 9 = 0.01 (T sa 1580 K), which yielded, respectively, the duration 
times, £ rca i = 5.35 seconds, 3.67 x 10 9 seconds, 1.81 x 10 21 seconds (pa 5.73 x 10 13 years), and 
2.52 x 10 79 seconds (pa 8.0 x 10 71 years). One notices tremendous duration times for the lowest 
temperatures. 

It is important to remind that the above results have not taken into consideration the de- 
generacy of the energy spectrum, which should contribute to decrease the duration of the quasi- 
equilibrium states (since the degeneracy increases the number of possible transitions taking 
place in the Monte Carlo process). However, the enormous estimates found above, for low 
temperatures, should not be altered significantly due to the degeneracy of the energy spectrum. 

B. Hydrogen Atoms in a Photon Bath 

It is well known that long-living nonionized hydrogen atoms exist in very low concentrations 
and at very low temperatures (typically 3 K), in intergalactic media. These atoms are in direct 
contact with photons, in such a way that transitions between states occur through photon 
emission and absorption. 

In what follows, we will present a crude estimate of a lower bound for the duration (in real 
time) of the above-mentioned quasi-equilibrium state for a hydrogen atom in a medium such as 
the intergalactic media. For that we will consider: 

(i) A nonionized hydrogen atom in contact with a photon bath at a temperature T; 

(ii) The temperature T sufficiently low for the atom to be found initially in its ground state; 

(iii) The velocity of the atom negligible with respect to the velocity of light. 

Let us define t\^ n {9) as the average time that the atom takes to absorb a sufficiently energetic 
photon [with an angular frequency u > A = (3R)/(4h)], in such a way as to perform a transition 
from the ground state to an excited state characterized by a quantum number n (n > 1), for 
a given value of 9. Obviously, t\^ n {9) represents a lower-bound estimate for the duration of 
the quasi-equilibrium state, since the atom may return afterwards to its ground state (this may 
occur with a large probability, since we are assuming low temperatures). Let us then consider 
the average number of photons per unit volume (including both directions of polarization), with 
angular frequency to > A, 
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«u;>a(0) = 



8vr 



uj 2 duj 



(2vrc) 3 A exp[(hu)/(R9)] - 1 ' 



(15) 



where we have used the variable 6 = ksT/R. For the temperature range of interest, one has 
that exp[(Hu)) / (R9)] » 1, in such a way that the above integral may be calculated easily, 



MS ~ 8 ^ f~ 3 

n ^ AW= (2^)3 exp [le 



2[?V + 2A^V + A^ 



h 



h 



(16) 



Using the result above one may estimate the number of photons with energy huj > 3R/4 in the 
volume of the hydrogen atom (to be considered herein as 47ra 3 /3, where a represents Bohr's 
radius). The maximum time that these photons spend within the volume of the atom is given 
by 2a/ c. Therefore, the average number of such photons in the volume of the hydrogen atom, 
per unit time, may be written as 



N^> A (0) = n w > A (0) (^) 



-na 



2a 2 

3vr C 2 6XP V 40 



R0 



+ 2A 



R0 



+ A ; 



R0 



(17) 



From this result one calculates the average time for the hydrogen atom to absorb a photon 
with sufficient energy to perform the transition 1 — > n, ti^ n (0)N UJ >\(0) = 1, i.e., ii_> n (#) = 
[-Wu>>a(0)] _1 - It is important to notice that, similarly to what happens for the duration of the 
quasi-equilibrium state - measured previously in Monte Carlo steps [cf. Eq. (5)] - the lower 
bound ii_> n (0) (in real time units) also follows an Arrenhius law, 



*i- n (0) ~ exp[3/(40)] = exp[3R/(4k B T)]. 



(18) 



However, as expected, the factor multiplying (1/6) in Eq. (18) is smaller than the one found in 
Eq. (5). One should notice that there may be alternative ways to obtain N u >\(9) through the 
knowledge of ?V,>a(0), as done in Eq. (17), e.g., by introducing a different time dependence in 
Eq. (17). However, the most important behavior, i.e., the Arrenhius law of Eq. (18), remains 
unchanged by using a different calculation for N u >a(0). It is important to mention that, in the 
present example there was no need for a connection between real time and Monte Carlo steps, 
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as done in the previous case. Even though, the low-temperature exponential-growth behavior 
appeared again; this result supports the hypothesis carried in the previous case for such a 
connection. 

Lets us now consider two typical examples for the lower bound ti^, n (0), namely 6 = 0.01 
(T ~ 1580 K) and 9 = 0.0001 (T ~ 15.80 K), which lead, respectively, to the colossal times (even 
when compared with the age of the universe of about 1.5 x 10 10 years), ti^ n (0.01) = 3.52 x 10 16 
years and ti^ n (0.0001) = 1.57 x 10 3243 years. Since the duration of the quasi-equilibrium state 
(in real time units) should be much larger than ti^ n (8), it becomes evident the treatment of the 
quasi-equilibrium state considered herein, as an effective equilibrium state for low temperatures. 

6. Conclusion 

We have shown that the hydrogen atom may live in a quasi-equilibrium state, at low temper- 
atures - in comparison with its corresponding ionization temperature - whose duration increases 
exponentially as the temperature decreases. By considering such a quasi-equilibrium state as an 
effective equilibrium state, we have calculated, for the first time (to our knowledge), thermody- 
namic properties within the BG statistical mechanics. For that, we have proposed a modified 
formalism (whose results are very close to those obtained through numerical simulations using 
the BG weight factor in the quasi-equilibrium state, at low temperatures), characterized by a 
regularized partition function. It should be stressed that such an approximation is supposed to 
be valid up to temperatures corresponding to about 6% of the ionization energy. Since the ioniza- 
tion energy E* of the hydrogen atom is extremely high (E* w 13.61 eV, i.e., E* /ks = 1.579 x 10 5 
K), our approximation should work well up to temperatures 10 4 K. It is important to mention 
the broad interest of the above analysis, which applies to many atoms, molecules, composite 
particles, and other similar systems, characterized by: (i) upper bounds, preceded by a quasi- 
continuum of energy levels, in their energy spectra; (ii) large gaps separating the ground and 
first-excited states. Obviously, experimental investigations are highly desirable in order to test 
the validity of these results. 
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